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We present a detailed analysis of the NA49 experimental particle yield results, and discuss the 
physical properties of the particle source. We explain in depth how our analysis differs from the 
work of other groups, what advance this implies in terms of our understanding, and what new 
physics about the deconfined particle source this allows us to recognize. We answer several fre- 
quently asked questions, presenting a transcript of a discussion regarding our data analysis. We 
show that the final NA49 data at 40, 80, 158 AGeV lead to a remarkably constant extensive ther- 
mal chemical-freeze-out properties of the fireball. We discuss briefly the importance of thermal 
hadronization pressure. 
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1. Deconfinement and Hadron Production 



Enhanced production of strange hadrons, and of (strange) antibaryons is a signature of quark- 
gluon plasma (QGP) formation in relativistic heavy ion collisions We illustrate on left in figure 
[T] the two step mechanism of hadron production from quark-gluon matter: upon deconfinement 
thermal glue emerges from parton matter, thermal gluon fusion reactions produce quark pairs in 
mass range m, <3T Subsequently, at a later time and thus at a much lower ambient tempera- 
ture, quarks merge into final state hadrons ^ ^. On the right in figure H we show the relative yield 
as function of p± of baryon to meson abundance which are quite different from what is found in pp 
reactions [^]. This confirms that the relativistic heavy ion hadron production differs from the string 
breaking mechanism of hadron production This result found also at top SPS energy range 




Figure 1: On left, qualitative illustration of (strange) hadron production in QGP: gluon fusion reactions 
populate flavor yield, subsequently at phase breakup, quarks merge into final state hadrons. On right, exam- 
ple of baryon to meson ratio experimental results A/Ks and p/ n in Au-Au compared to pp collisions as a 
function of Pi_. 



Of particular relevance in this work is the parameter controlling the absolute yield of quark 
pairs '^^^ for m, d and yP*^^ for s (and y^^^ ^ further discussed here). In an 

equilibrating system all y^'^^ — > 1 as function of time, one can argue that at RHIC strange flavor 
in QGP is nearly equilibrated. However, in an experiment, we observe hadrons, these yields are 
subject to different parameters operating in the hadron phase (HP) with the initial values 
If hadrons do not scatter much after formation, as is the case in the single freeze-out model |^, 
there is little if any change of yj^^itf). 

Remarkably, these two particle yield parameters, y^^{tf) and 7^^(?/), also control the relative 
yield of baryons to mesons shown on right in figure 



baryons 7™^^ f 7™ 



..HP I la ' 



mesons y™^ '"^^ ^ 



where yq shifts the yield of strange vs non-strange hadrons, thus for example the relative yields 
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of hadrons obey: 

A{ads) K+{us) 7™ <p a{sss) 72HP 

— oc — — oc — — oc — oc 

p{uud) 7^ n+{ud) 7™ h A{sud) 72HP 

The above observation, when combined with the resuhs shown on right in figure [l], imply that 
the baryon to meson yield and thus y^^{tf) is arising from microscopic dynamics of hadronization. 
This establishes the necessity to include the occupancy parameters in order to describe the yields 
of hadrons, since this is the parameter which allows for a reaction dependent relative yield of 
mesons and baryons. Conversely, a study of particle yields with a fixed value = 1 presumes 
certain relative yield of baryons to mesons and this is over-constraining a model of hadronization. 
Therefore, in such an approach, in order to describe the data precisely, the control of the relative 
number of produced mesons and baryon is achieved by manipulating the hadron spectrum [|l^] — 
including additional meson and baryon resonances, it is possible to fine tune the final relative yield 
of meson and baryons. Naturally, such a fine tuning of the hadron mass spectrum in fact confirms 
the failure of the 'equilibrium' model [11], and confirms the necessity to use the non-equilibrium 
occupancy parameters, if not for any fundamental reason, than in order to fudge the errors in the 
mass spectrum of hadrons. 

The greatest experimental sensitivity of the baryon to meson ratio is present when the respec- 
tive yields are comparable, such that the particle yields are a sensitive measure of this ratio. This 
is the case for the intermediate energy range at SPS, e.g., in the 80 AGeV Pb collisions with lab- 
oratory fixed target. Our choice to consider the SPS energy range in this presentation is in part 
dictated by this observation, and in part by the fact that we expect a change in reaction mechanism 
at an energy range below the RHIC 'high' energy runs. 

Other statistical hadronization (SH) parameters we derive from the data are the source volume 
V and the temperature T , at which particles stop changing in yield (chemical freeze-out). More- 
over, we obtain chemical potentials }Xb = 3/i^ = ?>T\nXq, /I5 = T\n{Xq/Xs), related to conserved 
quantum numbers: baryon number and strangeness, respectively. We also obtain A/3 which allows 
to conserve the net charge and expresses the asymmetry in the 3-rd component of the isospin. Es- 
pecially for low energy reactions where the particle yield is relatively low this parameter differs 
significantly from unity. 

The particle yields measure the values of SH parameters assuming that all particles stop chang- 
ing their abundance nearly at the same point in time. Here, the yield of stable hadronic particles 
(ignoring weak decays) includes hadronic resonance decay chain. Note that even if hadronic res- 
onances continue to react and their number changes in time after QGP breakup, the final number 
of stable particles is unchanged. Thus, the chemical freeze-out values of SH parameters, and the 
associated physical properties of the system at hadronization can be related to the QGP breakup 
condition, if these hadrons are free-streaming out of the deconfined fireball. 

Since the phase space density is in general different in the two phases, in order to preserve 
entropy (the valance quark pair number) across the phase boundary there must be a jump in the 
phase space occupancy parameters 7^ — this effect replaces the increase in volume in a slow 
re-equilibration with mixed phase which accommodates transformation of a QGP entropy dense 
phase into HP dilute phase. Similarly, there is a jump in strangeness occupancy since QGP is a 
more strangeness dense phase than is HP. y™{tf) > 7^'^^(?/), i = q, again due to rather large 
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Table 1: AGS (on left) and SPS energy range particle multiplicity data sets used in fits (see text). In bottom 
of table, we show the fitted statistical parameters and the corresponding chemical potentials. 
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0.017±0.005 


0.16±0.05 


0.39±0.06 


0.68±0.1 


1.82±0.36 


3.68±0.55 






1.5±0.3 


2.42±0.48 
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3.8±0.87 


4.5±0.20 
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0.14±0.07 




81±4 




3649±331 


4775±261 


2229±340 


1595±383 


2135±235 


3055±454 


T [MeV] 


153.5±0.8 


151.7±2.8 


123.8±3 


130.9±4.4 


135.2±0.01 


136.0±0.01 




5.21±0.07 


3.53±0.09 


2.86±0.09 


2.42±0.09 


1.98±0.07 


1.744±0.02 




1.565* 


1.39±0.05 


1.45±0.05 


1.34±0.06 


1.25±0.18 


1.155±0.03 
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0.49±0.03 
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0.877±0.05 
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276 
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180 


141 
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77.6 


62.0 


56.0 



hadron masses. In order to relate the initial values yj^^itf), i = q, s, at freeze-out (subscript-/), to 
the QGP fireball source value Y^^^itf), i = q, s, we match across phase boundary by considering 
the continuity of entropy and strangeness. In principle, only the occupancies are discontinuous 
at QGP hadronization, the other SH parameters are smooth. For this reason, we attached phase- 
indices 'QGP' and 'HP' only to occupancy parameters, and when these are omitted we always 
mean implicitly the confined hadron phase HP. 



2. Hadron Yields: Data and Fits 



A complete analysis of experimental hadron yield results requires a significant book-keeping 
and fitting effort in order to allow for resonances, particle widths, full decay trees, isospin multiplet 
sub-states. A program SHARE (Statistical HAdronization with REsonances) suitable to perform 
this data analysis is available for public use ^ 11]. This program implements the PDG [IM] 



confirmed (4-star) set of particles and resonances, and we use [ ]18[ ] the recent determination of a- 
meson mass ^] (m^ = 484, Fa/l = 255 MeV). The data set we use in table |I| is the latest NA49- 



2008 results [20], together with, as reference, AGS data set at their highest energy, see |21, p^ ]. 
A star at indicates that we fixed the value by strangeness conservation, which we do at AGS 
considering a small data set. 
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Table 2: Predictions: AGS/SPS particle yields including NA49 2008 data, all SHARE 2.2. 
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0.43 


0.69 
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Aside of directly measured yields, we fit the baryon content ('measured' in terms of event 
centrality choice), the charge per baryon ('independently' measured to be the proton content in 
nuclei) and the strangeness balance, all three shown in separate top data section of the table. Since 
'strangeness content' s + s is large, we choose to consider 5s = {s — s)/{s + s). The errors we 
present for these entries are our estimates of how the error in measurement propagates into the 
statistical parameters. 

A complete set of particle yields corresponding to the best parameters obtained at different 
reaction energies is presented in table ^. Below the stable particle yields, we show, in bottom 
section, the resonance yields which are not the measurable yields, but the initial values required 
in the study of a further resonance evolution |22]. The stable hadron yields agree well with the 
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Figure 2: On left: ;i^^/dof (top) and the associated significance level Pc[%] (bottom) as function of fixed 7^, 
for the AGS/SPS energy range. On right: the statistical parameters resulting from a fit at fixed 7^, as function 
of reaction energy. 



experimental energy dependence of the yield data well, this is, in particular, also the case for the 
-yield. However, we expect a growing yield of A, which experimentally saturates at 80 AGeV 
and then taking the center of measurements point, perhaps even decreases at 158 AGeV compared 
to80AGeV. 

We beheve that imposing strangeness conservation can over constrain the fit, a step we only 
take if the data sample is small so that there is no sensitivity to the 5s = constraint, as is the case 
at AGS. The values we find for ds at SPS are all negative, see values stated in top of table |. The 
persistence of 5s < 0, if not result of some experimental NA49 bias, should be taken as indirect 
evidence for some new physics. We are looking forward to RHIC low energy run to resolve this 
important question. Note that the particle yields presented in table ^ are evaluated with the value 
of 5s as found in fit to be the best value. All yields shown in this table are prior to weak decays 
after hadronization. 

We make, here, an effort to determine if the SH model we use makes good sense and if the 
parameters we employ are necessary. For this purpose, we show on left in figure ^ in the top 
frame the profile of ;f^/dof as function of 7^, and in bottom frame, the profile of the confidence 
level Pc(Z^)dof). We note that there is significant variation as function of 7^ and that the best fit 
is not associated with the value 7^ = 1. This shows that this parameter not used by many other 
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groups [11, 15, 16, 17| is neither redundant nor can it assume a fixed value, tacitly set to jq = 1, 
associated with thermal equilibrium ratio of baryons to mesons, (note that in the following, we 
always refer to and thus we sometimes omit as in this paragraph the superscript). 

We see, on left in figure that our fit results sepaiate into two groups. The AGS data point 
and the lowest energy SPS data point have a clear preference for < 1. This suggests that at 
these two reaction energies (11.6 and 20 A GeV) the yield of hadrons and in particular of baryons 
are suppressed, where the benchmark production level is a chemically equilibrated hadron phase at 
the temperature T ~ 152 ± 2 MeV. For the 30, 40, 80 and 156 A GeV SPS reactions, we see that 
our fit results favor enhanced yields of hadrons and in particular that of baryons compared to the 
benchmark which is a rather low chemical freeze-out T < 140 MeV. 

Sometimes it seems to us that other groups resist the use of 7^^ since they did not incorporate 
this important quantity consistently in their fit programs. Claims that SPS fits which include 7^ are 
unstable were not true for SPS energy range and are certainly false after the final complete data set 
is now available. Aside of 7^, there is another important nuance between our and 'their' fit. We do 
not see a value of A/3 published by several other groups, thus we are not sure that this important 
quantity is employed in the fits, along with the requirement that net electrical charge of the system 
is that of the participating protons. 

On right in figure ^ we show the fit results for statistical parameters, shown for each y^^NN. 
We note that there is a rapid shift in behavior as we go across ^/sun = 7 GeV (vertical dotted 
line) of r, Yq, and Jq/js, and of the volume V seen in the table |l|. This contrasts with the much 
smoother behavior of the baryo chemical Hb and strangeness /I5 potentials (bottom frame on right 
in figure ^ showing that across the entire reaction energy domain the chemical conditions change 
smoothly. On the other hand, there is considerable shift in chemical potentials /x, of individual 
baryons and anti baryons at the boundary ^/^nn = 7 GeV. Individual are evaluated following the 
quark content, for example: 

Ms(..s,) =2rin(7, + A,) + rin(7, + A,)), Msc^,-) = 2rin(7, - A,) + rin(7, - A,)). 

In figure ^ on left, we see that considerable discontinuity arises for all baryon chemical potentials 
between these two domains. The energy needed to add a baryon is thus smooth, while the energy 
to add any individual particle carrying baryon number is not. 

In figure ^ on right, we compare our fit of K+/ 71+ with data as function of ^/5nn- While our 
non-equilibrium model works very well and the peak of the horn we see in the data appears to be 
another consequence of the change in hadronization behavior near y^^NN = 7 GeV, the chemical 
equilibrium model [jlj] (dotted, blue) does not have a priory this capability. We believe that in 
order to explain the NA49 results one needs: 

a) the strangeness yield parameter 

b) and ability to choose baryon to meson ratio independent of temperature, that is 7^ 1 . 

The above parameters could be de-facto achieved by a fine tuned spectrum of hadron resonances 



[10]. We note that the decreasing value of K+/7r+ beyond the peak has recently been confirmed by 



the trial preparatory low energy run at RHIC [23] 



In closing this discussion, we note that an analysis of high energy RHIC data lacks the impor- 
tant constraint to a fixed baryon content, since unlike at SPS the RHIC results are at central rapidity 
and the baryon content is one of variables which are outcome of the fit. The absence of baryon yield 



7 



Particle Production and Deconfinement Threshold 



Johann Rafelski 



500 



S 

'— ' -500 



-1000 



-1500 

























* ... 






■ Mn 








•-- 


• - -. 


- «. 

- • 


— • 

-- • 


• Ma 








0-- 


--0- --- 


-: ::S - 


-----0 
b 


* Mz 












;::;:8;---- 




• Mn 
















o /in 




— 












- 






p 










° MX 
















° Ms 





































10 



0.25 



0.2 



0.15 



- + 


, 1 , 1 










- ^ I 










, , , 1 





[GeV] 



10 

^ [GeV] 



Figure 3: On left, the chemical potentials of all baryons and antibaryons in HP. On right, our model 
solid line (blue) tracks K+/;r+ yield (data) as function of reaction energy. Dotted (red) line: equilibrium 
model 

constraint reduces the stability of the many-parameter fit and a different approach must be taken 
which we will discuss elsewhere. 



3. Discussion of data analysis 

Q: As I see you do not constrain net strangeness to zero — I think that this maybe the reason your 
effort in this field is not taken as seriously as it should! 

A: Indeed, we have gained better insight in past 15 years regarding strangeness conservation, as 



compared to our older work [24]. At the time we argued that the slowness of the weak strangeness 
decay assures that in hadronic interactions net strangeness should be zero: As = s — s = 0. Today, 
we allow 5s to behave like a measurement, that is, we introduce it with an error 5s = 0±0.05. The 
reasoning is as follows 

a) This is a test of the hypothesis that weak decays remains weak in QGP phase, and the net 
strangeness remains conserved; 

b) NA49 did measure many, but not all particles carrying strangeness, e.g., has not been mea- 
sured, this yield is uncertain and thus 5s = cannot be tested to better than about 5%; 

c) Summing all measured and unmeasured hadrons in strangeness 'conservation' condition, = 
Li K ~ H ~^ 0» combines independent measurement errors and thus, even if NA49 had measured 
all strangeness carrying hadrons, there would be a residual statistical error present in 5s — another 
way to understand this is to note that we cannot confirm that weak decays in QGP remain weak to 
better than the progressing error of individual contributing measurements; 

d) Some strangeness could escape detection in unknown 'particles', for example being bound in 
(nearly) Mt/^-quark-symmetric semi-stable strangelett (a small drop of quark matter), this leads to 
5s <0 — which is what we find consistently in our fits. 

Q: In what sense is Pc (see figure^ left bottom frame) conveying confidence, and in 'what thing', 
in what way does this quantity deserve this name, and, I would think it is better to suppress this 
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part of the figure as it contains redundant information! 



A: You are familiar with the Section 32.3 in PDG [ ]14| ] (PDG: particle data group biannual issue 
of 'Review of particle physics'), which provides further discussion of the elusive meaning of Pq, 
and figure 32.2 shows which fixed values of Pq arise for given values x^/Aoi and dof. The SHARE 
program evaluates Pq for each fit. Abbreviating this section in PDG in a few words: When fitting 
a statistical data set to a model with a few parameters the value of Pc(;^^,dof) predicts (in sense 
of Bayesian inferment) the likelihood of repeat experiment to produce the same outcome for the 
model parameters. Thus, Pq expresses confidence in the validity of the model. 

Effectively, Pc(Z^)dof) also expresses confidence in the data, provided that we believe in the 
model: when we intentionally alter an experimental data point by 2 s.d., our data fit remains stable 
in the sense that we find nearly the same model parameters, but Pq becomes much smaller, and the 
one data point which is not well fitted is the one we altered. Note that ;^^/dof is not as sensitive 
to this consideration, since the effect of one 'wrong' measurement is diluted by the magnitude of 
dof. The function Pc(Z^)dof) turns insensitive in the limit dof oo, but we are far from this. In 
fact, in our study of the fits to the date, we found that the value of Pq is much more interesting 
when comparing for example NA49 data at different energies since the number of available data 
points varies and thus dof varies along with x^. Thus, if we were to suppress a result it would be 
the profile of ;^^/dof, but that cannot be done since people want to see the value of ;^^/dof. 



Q: I note that in your recent work [25] all 6 best values Pq are larger than 0.8 which is very unlikely 
to happen for 6 independent measurements. 

A: Yes, you will further note that in this paper all NA49 fits carried out with the final 2008 data 
set converge to a common best value Pq = 70%. This value has higher Pq than one could expect 
on statistical grounds (50%), but we note that we have treated systematic errors as if these were 
statistical errors, and thus effectively made errors too large. In one case, we also increased the error 
bar of NA49, beyond the statistical and systematic error, namely we doubled error for A-yield at 
80 AGeV, and thus it is of same relative magnitude as the error at 158 AGeV. We did this since our 
fit and Pc indicated that this one NA49 value was out of systematics of other particles measured at 
this energy. We fit this measurement but with an error which is modified as stated above. 

Q: How did the SPS heavy ion energy scan research program come to be? 

A: In early 1980's there were several studies of the mechanism governing the development of de- 
confinement as function of reaction energy. Looking at this work one finds energy estimates for 
transition to QGP within a range well beyond the SPS reach. Only a more persistent search will 
unearth a minority view, for example in the opening of conclusions of Ref.[26] we read: "The 
formation of a baryon-rich quark-gluon plasma appears to be an important reaction channel in col- 
lisions of heavy nuclei in the energy region of 2.5-5 GeV per nucleon in the center of mass frame of 
reference". We recall that the K+/ 7r+ peak is at 2 x 3.5 GeV, right in the middle of this prediction, 
obtained considering the growth of 'QGP seeds', small deconfined drops of matter. Strange particle 
production signature of this transformation was the favorite experimental approach [27|. However, 
only when it was proposed that an easily accessible observable K+/7r+ could suffice [28], the re- 
search program to scan the hadron production as function of reaction energy in the SPS domain 
took off. This development is consistent with the CERN announcement of February 7, 2000 that a 
"new state of matter" is produced in the top SPS energy central Pb-Pb collisions. 
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4. Physical Properties of the Source 

The SHARE program remains till this day a unique tool when it comes to the evaluation of the 
physical properties of the source of hadrons. Once a fit is achieved, this is a straightforward task, 
however, with a significant bookkeeping challenge: each hadronic particle, stable, or resonance, 
given the statistical parameters T, jj., contributes to the physical properties such as pressure P, 
entropy density S/V, energy density E/V, and baryons also to net baryon density b/V, while all 
strange hadrons with s-contents contribute to s/V. These properties are of considerable importance 
in discussion of the meaning of the fit to the data. 

We believe that the extensive properties such as P, or s are more reliable than the intensive 
statistical model parameters. The probability Pq is in fact indicating the reliability of finding in a 
redo of experiment and analysis in framework of the same physical model the same P, S, E, b, s. 
This is so since some elements of our model are not fully established and thus will evolve in 
time. For example, the hadron spectrum will be in the coming decade better understood and this 
may alter some statistical properties. However, the occupancy parameters jj compensate to a large 
extent just this model dependence, and thus a refinement of hadron mass spectrum should result 
in some change of ji, but, as we hope, little changed evaluation of the physical properties of the 
fireball. Said differently, for example, the observed pion yield is a good measure of the entropy 
produced with little model dependence. 

A sample of our findings for the physical properties is presented, in figure ^, as ratio of differ- 
ent extensive quantities: on left, we see s/b, s/S, and £'/5'; on right top frame, the pressure P, and 
bottom right, the thermal energy E per primary hadron hp. The results we find are the symbols, the 
lines guide the eye. All solid (blue) symbols and lines are full non-equilibrium model. On right 
the open, dashed (red) triangles and lines are semi-equilibrium model with = 1. Since y^^ 
is artificially fixed, the pressure is not the same in both models, and, there is a great reduction of 
confidence level Pq in this result, as is seen in figure ^ left bottom frame: compare highest value 
of Pc with that arising at a fixed 7™ = 1. Said differently, the values we present on right in figure ^ 
for 7^^ = 1 (red, dashed) and 7^^ fitted (solid, blue) are different since these are in principle two 
different models to be distinguished by the quality of the fit. 

Perhaps the most interesting finding is seen in figure^ on right in the top frame: the hadroniza- 
tion pressure is practically constant, near about 82 MeV/fm^ (non equilibrium model, solid line), at 
highest three SPS energies. The absence of fluctuations in this result is most remarkable, a smooth 
line is connecting 40, 80 and 158 AGeV result. The preliminary NA49 data gave less 'constant' 



result [g5|,|18|]. 

The baryon content of the source fireball is fixed by the choice of centrality, thus the rise 
of s/b (top left) indicates that strangeness production increases rapidly, and smoothly, with re- 
action energy. Note that strangeness content increases six fold between 10 and 158 AGeV reac- 
tion energy, and by a factor 2.5 between 30 and 158 AGeV. On the other hand, the value range 
s/S G (0.019,0.022) combined with the growth of s/b shows that the fireball experiences, as func- 
tion of reaction energy, a growth in both entropy and strangeness at a comparable rate. 

The high value of 5/5 is established much earlier in time, in a much more dense phase than 
the chemical freeze-out. The behavior of s/S suggests that at 40 AGeV (and perhaps event at 30 
AGeV), in an early stage of the fireball when s and S is made, the same quark-gluon degrees of 
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[GcV] Vs^ [GeV] 

Figure 4: Physical properties of the hadronizing fireball. On left 4s /b, s/S, E/s and in right pressure P 
and thermal energy per primary hadron E /hp. The red Unes on right correspond to — 1 and have a much 
lower confidence level, Pq, see figure 2. 



freedom are already active. Consideration of kinetic theory [ p9[ ] suggest that the growth of 5'/5' is 
due to a gradual as function of reaction energy approach to chemical equilibrium, first by gluons, 
than by light quarks, and later by strange quarks. We refer to a more thorough discussion of the 



importance of observable 5/5 presented elsewhere [30, 31]. 



The cost in global thermal energy E to make a strange quark pair yield s, including those bound 
as pairs in T] , T]' and is seen on left, in the bottom frame of figure ^. This value is as low as 6 GeV 
at highest SPS energy. Here, E is the thermal energy content of the hadron system and s is as before 
the absolute yield of strange quark pairs. A strange quark pair at hadronization has an energy near 
to 0.9 GeV. The remainder of the energy content per this pair, in a nearly equilibrated deconfined 
source of hadrons, is contained in the thermal light quarks and gluons. This argument accounts 
well for the remaining 5 GeV. The greater energy required to produce strangeness pair at the lower 
reaction energies is consistent with the hypothesis that, prior to hadronization, strangeness in the 
fireball is far below chemical equilibrium, as is seen in the small s/b ratio. 

Finally, the bottom frame on right considers the hypothesis that the chemical freeze-out con- 
dition is related to energy required for making a primary hadron [|l^]. We see that within the 
semi-equilibrium model (red, dashed line), the value is indeed, within a reasonable margin, oscil- 
lating near E/hp ^ I GeV. However, the high confidence result (blue, solid line) within the non 
equilibrium model demonstrates a strong variation of E/ hp which drops from above 1.3 GeV to 
0.75 GeV, and thus, E/hp is not a good criteria to define freeze-out of hadrons. Moreover, this 
variable has a priory no fundamental meaning, and this quantity depends strongly on the spectrum 
of hadrons used, equivalently, the value of 7^. 
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5. Discussion and Conclusions 

For low reaction energies (AGS 10.6 AGeV and SPS 20 AGeV), we find hadrons originate from 
a relatively large {V > 3500fm^), dilute (7^ < 0.5), yet relatively hot system {T > 150 MeV). The 
strangeness abundance is low and fast growing. At higher reaction energies (SPS 30, 40, 80 and 
158 AGeV), the volume per hadron is less than half as large, and the chemical freeze-out proceeds 
at about 20% lower T. We know from the experimental results that there is a rapid expansion of 
the dense matter. We conclude that in the buildup of the matter flow there was significant cooling 
of the fireball matter. Particle spectra confirm that there is fast flow in this SPS energy range. Our 
hypothesis, which is consistent with all results described, is that we are, already at SPS, seeing 
sudden hadronization of a supercooled deconfined phase [^] at 40, 80 and 158 AGeV reaction 
energy, and to some lesser extend at 30 AGeV. This hypothesis allows to qualitatively explain the 
high > 1 value, and is consistent with physical properties of the system we observe implicitly 
by measuring the hadronization parameters. 

There are several possible explanations of the AGS 10.6 AGeV and SPS 20 AGeV low energy 
behavior: 

a) the most straightforward alternative is the usual confined, hadron phase in which antibaryon an- 
nihilation, which reduces the baryon to meson ratio, takes place, which leads to chemically under 
saturated 7^ < 1 ; 

b) given the large baryo-chemical potential (see bottom of tablejl]) the freeze-out is in a domain 
in which a complicated QCD phase structure is expected and the physical properties we find are 



consistent with breakup of a valance quark deconfined matter [|18|]. Such a phase has relatively 
massive constituents (niq ~ 330, ~ 500 MeV) chiral symmetry is strongly broken, even though 
deconfinement prevails — only for /x, there is empirical and lattice evidence that chiral sym- 
metry and deconfinement are restored at the same condition. The matter flow of a heavy constituent 
matter is expected to be slow, and thus, hadronization mechanism differs profoundly from what one 
sees at higher energy; 

c) Considering relatively high Pq near 7^ = 1, for 10.6 AGeV we have 45% and for 20 AGeV 
we have 60%, we could simply discount our results at low energies forcing chemical equilibrium 
Yq = Ifor light quarks, claiming that the (small) best fit values of jq are an artifact of experimental 
data points. To some extent this remark can be made also for 30 AGeV data, but cannot be made 
for 40, 80 and 158 AGeV. 

Perhaps the most intriguing result of this analysis is the smoothness, and even near constancy, 
of physical properties of the fireball at chemical freeze-out condition for the top three SPS energies 
40, 80 and 158 AGeV. Of particular physical interest is the value of hadronization pressure P ~ 82 
MeV/fm^, which is nearly constant. In a phase transformation from quarks to hadrons, the pressure 
of quarks is transferred into the pressure of color-neutral hadrons, which can escape from the 
deconfined fireball. Since the flow pressure of quarks transfers smoothly into that of hadrons, we 
conclude that the thermal pressure of produced hadrons, /" ~ 82 MeV/fm^, provides a first estimate 
of the pressure of the vacuum which keeps color charged quarks inside the fireball up to the point 
of sudden fireball break-up. 

To conclude, we find analyzing the final NA49 SPS data that the most central Pb-Pb collisions 
at 40, 80 and 158 AGeV, corresponding to ^^nn = 8.76, 12.32, and 17.27 GeV, are showing very 
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similar behavior at hadronization, where the physical and statistical properties of the source can 
be interpreted in terms of a sudden QGP fireball decomposition. The SPS reaction at 30 AGeV, 
corresponding to ^.vnn = 7.61 GeV, is an intermediate case, which in many aspects is similar to 
the higher energies, but the hadronization pressure evaluated in the ensemble of reaction events is 
significantly reduced. The reactions at AGS 10.6 and SPS 20 AGeV, corresponding to y^Snn = 4.84 
and 6.26 GeV, are in a different class which cannot be yet fully categorized, and require more 
precise and comprehensive experimental data in order to allow a more specific conclusion about 
the appUcable reaction mechanism. 
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